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QJ , Abstract 

The fermion determinant and the chiral anomaly of lattice Dirac 
K> , operator D on a finite lattice are investigated. The condition for D 

5h I to reproduce correct chiral anomaly at each site of a finite lattice for 

smooth background gauge fields is that D possesses exact zero modes 
satisfying the Atiyah-Singer index theorem. This is also the necessary 
condition for D to have correct fermion determinant ( ratio ) which plays 
the important role of incorporating dynamical fermions in the functional 
integral. We outline a scheme for dynamical fermion simulation of lattice 
QCD with topologically proper D. 



PACS numbers: 11.15.Ha, 11.30.Fs, 11.30.Rd 



1 Introduction 

The fermion determinant, det(D), of a Dirac operator D, is proportional to 
the exponentiation of the one-loop effective action which is the summation of 
any number of external sources interacting with one internal fermion loop. It 
is one of the most crucial quantities to be examined in any lattice fermion 
formulations. In numerical simulations of lattice QCD, its ratio between two 
successive gauge configurations determines whether the effects of dynamical 
quarks are properly incorporated, which are believed to be cruical to bring the 
ab initio predictions of lattice QCD to agree with the experimental data. An- 
other important quantity in QCD is the chiral anomaly which breaks the global 
chiral symmetry of massless QCD through the internal quark loop. In the case 
of the axial current coupling to two external photons, the axial anomaly |I[ 
provides a proper account of the decay rate of 7r° — ► 77. For the flavor sin- 
glet axial current coupling to two gluons, the presence of chiral anomaly but 
no associated Goldstone boson posed the Ua(X) problem which was resolved 
by 't Hooft 0, after taking into account of the topologically nontrivial gauge 
configurations, the instantons. This may also provide an explanation why the 
77 ( rf ) particle is much heavier than the 7r's ( K's ). Lattice QCD should 
be designed to provide nonperturbative and quantitive answers to all these 
and related problems. However, these goals could not be attained if the lat- 
tice Dirac fermion operator does not reproduce correct chiral anomaly or the 
fermion determinant ( ratio ) on a finite lattice without fine-tuning any param- 
eters. Furthermore, lattice QCD should realize the spontaneous breaking of 
chiral symmetry in massless QCD, which gives massless Goldstone bosons and 
thus leads to a correct and detailed description of the properties of the strong 
interactions at low energy ||. 

A basic requirement for lattice Dirac operator D is that, on a finite lattice 
with any prescribed^] smooth background gauge field which has a well-defined 
topological charge, D possesses exact zero modes satisfying the Atiyah-Singer 
index theorem ( i.e., D is topologically proper ) without fine tuning any pa- 
rameters. Here we would like to emphasize the importance of having exact 
zero modes on a finite lattice without fine tuning any parameters. If the zero 
modes can only be obtained by fine tuning parameters, then it is impossible to 
tune the parameters for every gauge configuration generated during a Monte 
Carlo simulation. 

For any D satisfying this requirement, the sum of the anomaly functionP] 
over all sites must be correct, and is equal to two times of the topological charge 
of the background gauge field. If it turns out that the anomaly function of D 
at some of the sites do not agree with the Chern-Pontryagin density, we can 



1 An example of prescribed background gauge field is given in the Appendix. 

2 Here we have adopted the terminology used in Chapter 22 of ref. H. 



perform the following topologically invariant transformation^ 

T{R): D' = T(R)[D] = D(1 + RD)- 1 (1) 

with some operator R such that D' is local, then the anomaly function of D' 
would be in good agreement with the Chern-Pontryagin density at each siteQ. 
It suffices to choose R to be a hermitian operator which is local in the position 
space and trivial in the Dirac space. The terminology topologically invariant, 
refers to the property that the zero modes and the index of D are invariant 
under the transformation ([Jj. Note that it is not necessary to perform any 
fine tunings of R since any local D' would give the correct chiral anomaly. 
Evidently, the set of transformations, {T(R)}, with the group multiplication 
defined by 

(T(R 1 ) o T{R 2 ))[D] = T(R 1 )[T(R 2 )[D}}, (2) 

form an abelian group with the group parameter space {R}, due to the follow- 
ing basic property, 

T(R 1 )oT(R 2 ) = T(R 1 + R 2 ). (3) 

In the topologically trivial sector, D~ l exists and ([]]) becomes 

D'- 1 = D- 1 + R. (4) 

Then there exists a unique 

# c = -^{zr\ T5 } T5 (5) 

such that 

D' 1 = D- 1 + R c (6) 

is chirally symmetric, i.e., 

15 d; x + D- X l5 = 0. (7) 

From the physical point of view, we must require the existence of a chirally 
symmetric D c = T(R C )[D] such that in the free fermion limit D c (p) — > i^f^p^ 
as p — > 0, otherwise, D is irrelevant to massless QCD. In general, we assume 
that there exists 

D^Dil + RcD)- 1 (8) 



3 The form of Eq. (fil) is similar to the general solution of Ginsparg- Wilson relation given 
in ref. H]. However, the transformation (]l]) can be used for any Dirac operators. 

4 Here we have assumed that the size of the finite lattice is large enough such that 
the finite size effects can be neglected, i.e., the size of the lattice is much larger than the 
localization length of D'. 



such that 

£c75 + 75^c = 0. (9) 

If D(D C ) is topologically proper, then D(D C ) must be free of species doubling, 
thus D c is non-local, according to the Nielson-Ninomiya no-go theorem ||. 
Furthermore, D c has singularities in topologically nontrivial background gauge 
fields [0]. Substituting ([8]) into (|9|), we obtain 

D l5 + l5 D = -2DR cl5 D. (10) 

This is the rejuvenated Ginsparg- Wilson relation M, In general, for any lattice 
Dirac operator D, if the chirally symmetric D c = T(R C )[D] exists, then D 
must satisfy the Ginsparg- Wilson relation. From this viewpoint, the Ginsparg- 
Wilson relation really does not specify the important attribute of a Dirac 
operator, namely its topological characteristics 0. Only for D is topologically 
proper ( e.g., the overlap-Dirac operator || ), D ( or its transform D' = D(l + 
RD)~ X ) can be used for lattice QCD, then the attractive features pointed out 
in refs. [|1(], [TlJ can be realized on a finite lattice. If D is topologically proper 
but nonlocal, the transformation ([!]) can be used to obtain a local D' such 
that the anomaly function at each site of a finite lattice could agree with the 
corresponding Chern-Pontryagin density for smooth background gauge fields. 
Moreover, the locality of D' also implies that its fermion determinant ratio ( 
with the zero modes omitted ) would agree with the corresponding value in 
continuum. 

2 The Anomaly Function 

In general, we consider lattice Dirac operator D which breaks the continuum 
chiral symmetry according to 

D l5 + l5 D = B (11) 

where B is a generic operator which is usually taken to be irrelevant ( i.e., 
vanish in the limit a — > ). However, we can construct the chirally symmetric 
action 

A s = $D a il) (12) 

using the chirally symmetric part of D, 

D a = hip-^D^). (13) 



Here the Dirac, flavor and internal symmetry indices are all suppressed. Then 
A s has the usual chiral symmetry, and the divergence of the associated Noether 
current^ is 

d»Jfa) = ^7 5 W)* + $D) xl ^ x - \i^B) x ^ x - ^ x (Bi/j) x (14) 

which satisfies the conservation law, 

]T^J*(x) = 0. (15) 

X 

As usual, the lattice is taken to be finite with periodic boundary conditions, 
and d^J^x) is defined by the backward difference of the axial current 

Vj(*) = £W(*Wj(*-A)] (16) 

such that it is parity even under the parity transformation, and the conserva- 
tion law Eq. ([T^) is also satisfied. 

If D does not possess exact zero modes in the background gauge field , then 
the fermionic average of d^J^x) can be evaluated as 

(0"J») = |/[#][#]^4(x)exp(-^^) (17) 

Z = /"[#][#] exp(-$£ty) (18) 

Using (|I4"D, we obtain 

d^JHx)) = ~ tr [(BD-^ix^) + {D- l B){x,x)] = A n {x) (19) 

where the trace runs over the Dirac, flavor and internal symmetry indices. 
The RHS of Eq. ( |T9| ) is identified to be the anomaly function of D. If the 
chirally symmetric limit of ([!]), i.e., D c = T(R C )[D] exists, then D satisfies the 
Ginsparg- Wilson relation flrO|) and the anomaly function becomes 

An{x) = -tr[ l5 (R c D)(x,x)+ l5 (DR c )( Xl x)} (20) 

However, for any D which does not possess exact zero modes in the background 
gauge field, the sum of the anomaly function ( [L9|) over all sites must vanish 
due to the conservation law, Eq. (|i~5|), 



EA»(ar) = 0. (21) 

X 

This implies that ii A n {x) is not zero identically for all x, then it must fluctuate 
from positive to negative values with respect to x. The latter case is exactly 

5 Here we assume that it is the flavor singlet of the flavor symmetry group. 



what happens to the anomaly function of the standard Wilson-Dirac fermion 
operator. 

On the other hand, if D possesses exact zero modes in topologically non- 
trivial background gauge fields, then D~ x is not well defined. In this case, one 
needs to introduce an infinitesimal mass m and evaluate ( |T7D with D replaced 
by D = D + mf[D] ( where f[D] is any functional of D, which has eigenvalue 
one for the exact zero modes of D ), and finally take the limit ( m — > ), i.e., 

yj%x)) = hmi|[#][#]^Jj(x)exp(-^^) (22) 

exp(-<0L>V) (23) 

Inserting fll4[) into (^), we obtain 

'd^Jlix)) = Um]-tT\(BD- 1 )(x,x) + (D- 1 B)(x,x) 

N + V_ 

+2^[^+(x)]t^(x) - 2 5>r(*)] t *r(*) (24) 



s=l t=l 



where 0+ and 0^ are normalized eigenfunctions of D with eigenvalues X s = 
Xt = and chiralities +1 and —1 respectively. The first term on the RHS of 
Eq. (El) is identified to be the anomaly function of D, 



An{x) = lim - tr 

m—>0 2 



(BD- 1 ) {x, x) + (D- l B) (x, x)] (25) 



Then summing Eq. fl2~4|) over all sites and using Eq. ([15|), we obtain 

iV_-iV + = i^X(x) (26) 

This is the index theorem for any lattice Dirac operator on a finite lattice, in 
any background gauge field. In the case D does not possess any exact zero 
modes, Eq. (^) reduces to Eq. (pTj) . Equation (|26|) implies that the sum 
of the anomaly function over all sites must be a well defined even integer for 
any background gauge field, however, it does not necessarily imply the Atiyah- 
Singer index theorem for smooth background gauge fields. This can be seen as 
follows. Since D does not have exact zero modes in a trivial gauge background, 
the index of D must be proportional to the topological charge Q of the smooth 
background gauge field. Now, if Q is an integer, then the proportional constant 
must be an integer, otherwise their product in general cannot be an integer. 
Denoting this integer multiplier by c[D], we have 

AL - iV+ = c[D] Q n f (27) 



where nj is number of fermion flavors. In particular, for n/ = Q = 1, 
c[D] = N- — N + . Here we have assumed that c[D] is constant for smooth 
background gauge fields. This is a reasonable assumption since c[D] is an in- 
trinsic characteristics of D. However, when the gauge field becomes rough, we 
expect that Eq. (|2"7| ) would break down. If one insists that Eq. fl27|) holds even 
for rough gauge configurations, then c[D] cannot be an integer constant due to 
the highly nonlinear effects of the gauge field. It is easy to deduce that, in gen- 
eral, c[D] is a rational number functional of D, which in turn depends on the 
gauge configuration, but it becomes an integer constant only for smooth gauge 
configurations. The topological characteristics of D, c[D], was first discussed 
in ref. |7j], and was investigated in ref. [12] for the Neuberger-Dirac operator. 



Equation (|2"7 ) constitutes the index theorem for any lattice Dirac operator on 
a finite lattice, in any background gauge field with integer topological charge. 
We can classify D according to its response in a smooth nontrivial background 
gauge field. If D does not possess any zero modes, then c[D] = 0, D is called 
topologically trivial. If c[D] = 1, ( ]2"7| ) becomes Atiyah-Singer index theorem, 
then D is called topologically proper. If c[D] is not equal to zero or one, then 
D is called topologically improper. 

Now we try to obtain a general expression for the anomaly function sat- 
isfying (|26|) and (|27|) . Consider the gauge configuration with constant field 
tensors, e.g., F^ 2 (x) = j^f-, Fu( x ) — "jrf 2 " an d other F's are zero, where qi 
and qi are integers. Then the topological charge of this configuration is 

Q = ^T2 E e ^- F^(x)F b Xa (x +fi + P) tr{t a t b } = n gi q 2 (28) 

O.Z7T -. 

which is an integer. Here t a and £& are generators of the internal symmetry 
group with the normalization tr{t a tb} = n 5 a b- If D is local, then it is rea- 
sonable to assume that A n (x) is constant for all x. From Eqs. (PB|)-(PH[), we 
obtain 



An(x) = ^ c[D] e, vXa F^(x)Fl(x + jx + u) ti{t a t b ] . (29) 

However, this implies that any local D would yield the correct anomaly func- 
tion on a finite lattice, up to the integer factor c[D]. In reality, this cannot 
be true in general. The crux of the problem is due to the fact that when D 
is topological trivial, i.e. c[D] = in the index theorem (27), the anomaly 



function A n (x) does not necessarily vanish identically at each site, but it can 
fluctuate from positive to negative values provided that its sum over all sites 
is zero. This is the complementary solution to the homogeneous equation 
J2 x A. n (x) = 0. Using Eq. (|l^), we can write the general form of the comple- 
mentary solution as 

A {h \x) = J2 [Gp{x) ~ G,(x - fi)} = d,G,(x) (30) 



where G^ix) is any local and gauge invariant function. Then the general 
solution of the anomaly function is the sum of the complementary solution 
and the particular solution 



An{x) = Y^ c[D]e ^ F ^ {x)F ^ {x + fi + i))tl{tatb} 

+ E [G>(aO - G„(a: - £)] (31) 

In general, G^x) depends on the background gauge field as well as the lattice 
Dirac operator D. If A^ (x) does not vanish in the presence of a constant gauge 
background, then A n (x) must be different from the continuum chiral anomaly, 
hence this D should be abandoned. Only when AS h '(x) vanishes identically for 
any smooth gauge background, D can reproduce the correct chiral anomaly 
for c[D] = 1 ( e.g., the Neuberger-Dirac operator ). On the other hand, 
the pathologies of the standard Wilson-Dirac operator are : c[D] = and 
A^(x)^0. 

Now we consider D with anomaly function A n (x) satisfying (|29|) in a con- 
stant gauge background. Then we introduce local fluctuations [ e.g., the sinu- 
soidal terms in Eqs. (|63|)-(|66|) ] to the background gauge field and/or increase 
its topological charge Q. Then c[D] in Eqs. (|29|) and fl27|) may remain con- 



stant provided that the local fluctuations of the gauge fields are not too violent 
and/or Q is not too large. When the gauge background becomes so rough that 
D is non-local, then Eq. (E9) is no longer valid, but the index theorem, Eq. 
(p7|) may still hold with the same c[D]. In this case, the anomaly function 
A n (x) can be expressed in terms of the general formula fl3"T| ) with the same 
c[D] but nonvanishing <9 M G M (x). If we keep on increasing the roughness of the 
gauge background, then D will undergo a topological phase transition, and 
c[D] in Eqs. (^7\ ) and (|3~ID will become another integer or even a fraction ( 
Explicit examples of topological phase transition will be demonstrated in Sec- 
tion 4 ). In general, Eq. ([H]) constitutes the anomaly function for any lattice 
Dirac operator Dona finite lattice, in any background gauge field. However, 
only when A- h '{x) = and c[D] = 1, the anomaly function can agree with the 
Chern-Pontryagin density. 

If the chirally symmetric limit of D, i.e., D c = T(R C )[D], exists, then 
D satisfies the Ginsparg- Wilson relation (|10D, and its anomaly function (p5|) 
becomes exactly in the same form of Eq. (|20|) . Substituting Eq. fl20D into the 
LHS of Eq. (|3T|), we obtain 

An{x) = -tr[ l5 {R c D){x,x)+j 5 (DR c )(x,x)} (32) 

' ' c[D] e, uXa F a (x)Fl(x + jl + v) tr{t a t b } 



lQ n 2 L J -vwv \lv\ 

+ J2lG,(x)-G ll (x-fi)} (33) 



which agrees with the results obtained in ref. |7|, [H|. For U(l) lattice gauge 
theory with "admissible" background gauge configurations, Eq. (^) agrees 
with the result derived by Liischer jTJ]]. For the Overlap-Dirac operator, in 
the weak coupling limit or the classical continuum limit, we have d^G^(x) = 0, 
then Eq. (^) agrees with the results derived by Kikukawa and Yamada flT5| , 



Fujikawa jTIJ, Adams (T^J and Suzuki |18| respectively 



Inserting Eq. (|25|) [ Eq. fl20|) ] into Eq. fl26|) , we obtain 

N- - N + = - l - £ tr [j 5 (R c D)(x } x) + l5 (DR c )(x, x)\ (34) 



which agrees with the index theorem for Ginsparg- Wilson fermion |19[ |20 . 
Note that (^) does not necessarily imply the Atiyah-Singer index theorem for 
smooth background gauge fields with integer topological charge. 

If the lattice Dirac operator D does not have exact zero modes in the pres- 
ence of topologically non-trivial background gauge fields, however, it manages 
to possess exact zero modes in a certain limit, say, for example, the Wilson- 
Dirac fermion operator || in the classical continuum limit ( i.e., zero lattice 
spacing with infinite number of sites in each direction ). Then on any fi- 
nite lattice, this lattice Dirac operator cannot reproduce the correct anomaly 
function for topologically non-trivial background gauge fields, since the sum 
of the anomaly function fllTf ) over the entire lattice must be zero due to the 
conservation law flHf), 



W tr{(BD- l )(x,x) + (D- 1 B)(x,x)} = £(«9^(x))=0 (35) 

* x x 

while the sum of Chern-Pontryagin density over all sites is proportional to the 
topological charge, 

E 7^2^ F ^)Fl{x + (i + 0) tr{t¥>} = 2 Q (36) 

which is nonzero for topologically nontrivial background gauge fields. As we 
increase the number of sites or decrease the lattice spacing, the anomaly func- 
tion does not approach the Chern-Pontryagin density since it must satisfy the 
constraint Eq. (p5|). It is evident that the limit where D might possess exact 
zero modes and correct anomaly function cannot be obtained by smoothly ex- 
trapolating measurements at finite lattice spacings towards a = since the total 
anomaly $35] ) must undergo a discontinuous jump from zero to a nonzero inte- 
ger at the point a = 0, if the background gauge field is topologically nontrivial. 
Only for the trivial background gauge field, D might have the chance to re- 
produce the correct anomaly by extrapolations towards the limit a = 0. If one 
attempts to bypass this difficulty by adding a mass term to D with the hope 
of obtaining the correct anomaly by tuning the mass and/or other parameters, 

8 



then one must encounter the unsurmountable problem of fine-tuning, since it 
is almost impossible to tune the parameters to some fixed values such that 
the Chern-Pontryagin density is reproduced at each site, and for all different 
gauge configurations. Therefore the chiral limit of these lattice Dirac operators 
( e.g., the standard Wilson-Dirac fermion operator ) not only is difficult to be 
realized in practice, but also in principle a goal impossibly to be attained. On 
the other hand, if the lattice Dirac operator ( e.g., the overlap-Dirac operator 
U ) can reproduce exact zero modes satisfying the Atiyah-Singer index the- 
orem on a finite lattice, then the total anomaly is bound to be correct, thus 
the anomaly function could agree with the Chern-Pontryagin density even at 
finite lattice spacings. 

3 Fermion Determinant 

Now let us turn to the fermion determinant which plays the important role 
of incorporating dynamical fermions in the functional integral. Formally, the 
fermion determinant can be represented as the product of all eigenvalues of D. 
In continuum, the fermion determinant must be regularized, and only the ratio 
of fermion determinants can be well defined and independent of regularization. 
In the topologically trivial sector, it is ususally normalized to one in the free 
fermion limit, i.e., det[-D(A)] <— det[Z)(y4)]/det[D(0)]. However, in the non- 
trivial sectors, D has exact zero modes, then the fermion determinant ratio is 
well-defined only for gauge configurations within the same topological sector 
or with the zero modes omitted. 

On a lattice, if D is topologically proper, then the fermion determinant in 
the functional integral also serves another important role to "decompose" con- 
figurations of link variables into distinct topological sectors. This can be seen 
as follows. On a lattice, any two gauge configurations are smoothly connected. 
In Monte Carlo simulations, the transition probability between any two gauge 
configurations in the quenched approximation is proportional to the exponen- 
tiation of the difference of their actions, which is finite. This is in contrast to 
the situation in the continuum, where the gauge configurations are decomposed 
into disconnected topological sectors, and the functional integral is performed 
in each topological sector individually. However, if D is topologically proper 
and its fermion determinant is incorporated in the functional integral, then the 
transition probability between any two gauge configurations in two different 
topological sectors is either zero or infinite, in the limit m — > ( Here m is an 
infinitesimal mass introduced to D ). So, if one cutoffs all transitions having 
a transition probability less than a certain lower bound, say, pi(m), or larger 
than a upper bound, say, p 2 (m), then only gauge configurations within the 
same topological sector can be accepted in the Monte Carlo updatings, thus 
the functional integral can be evaluated exactly in the same manner as in the 



continuum with disconnected topological sectors. Therefore Monte Carlo sim- 
ulations and measurements of observables are performed separately for each 
relevant topological sector, and the final result is averaged over all relevant 
topological sectors with weight e lud , where 6 denotes the vacuum angle and 
v the winding number. From this viewpoint, the deficiency of link variables 
being smoothly deformable to the identity element can be circumvented by 
incorporating dynamical fermions in the functional integral, provided that D 
is topologically proper. 

To be specific, consider lattice QCD with nf flavors of "massless" quarks, 
where the lattice Dirac operator D is topologically proper. Then the expec- 
tation value of a general observable 0(U, ip, iji) of the field varibales can be 
written as 

(O) = i J[dU][d^}[d^]0(U,ij,tP)exp(-A g -i>D(U)^) (37) 

where A g is Wilson's pure gauge action, and Z is the partition function 

Z = f[dU}det[D(U)}exp(-A g ) (38) 

Here the functional integral of gauge fields must include the contributions of 
instantons. In four dimensions, for any gauge group G which is a simple Lie 
group ( e.g., G = SU(3) ), the instantons are characterized by the elements of 
the homotopy group vr 3 (G) = n 3 (SU(2)) = Z, i.e., the winding number, 

v = q^2 £ ^x«F; u {x)Fl{x + (i + i>) (39) 

which satisfies 

n f v = AL - iV + (40) 

where N + (N^) denotes the number of zero modes of D of positive (negative) 
chirality. Suppose that the weight factor for each topological sector is w(is), 
then Eqs. ( |T7D and Q3"8|) are replaced by 

(O) = ^Y, w ^)Jl dU }[dnd^O(U^^)ex V (~A a -^D(U)i;) (41) 

Z = J2 W ( U ) f[dU]det[D(U)]exp(-Ag) . (42) 

If we require (O) to satisfy the cluster decomposition principle ||, then it is 
straightforward to deduce that 

w{v) = exp(iu6) (43) 
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where 9 is a free parameter, the so called vacuum angle pl| . Due to the 
presence of fermion determinant, the partition function fl4~2|) only receives the 
contribution from the trivial sector, 

Z = [[dU]( 0) det[D(U)] exp(-^ 9 ) (44) 

where the superscript (0) stands for the trivial sector. Without loss of gener- 
ality, we write the observable O as 

0(U, i>, i>) = 0{U)^ Xl i, yi ■ ■ ■ ^ x J yn . (45) 

Then, for the trivial sector, ([H]) can be evaluated as 

(0) = y[dUreM-A g )0(U)det(D) £ e^;- D^ yi ■ ■ ■ D^ (46) 

where e^ 1 '.'.'.^™ is the completely antisymmetric tensor which is equal to +1(— 1) 
if {zi ■ ■ ■ z n } is even ( odd ) permutation of {xi ■ ■ ■ x n }, and is zero otherwise. 



For topologically non-trivial sectors, (|4"TD can be evaluated by expanding the 



fermion fields in terms of the complete set of eigenfunctions of D with Grass- 
man coefficients, then the functional integral over the fermion fields can be 
expressed in terms of that over the Grassman coefficients. It is evident that 
for O given in Q4"5]), non-trivial sectors with winding number \v\ > n do not 
contribute to (|4i"|). In general, ( f4i~D can be written as 

(°) = ^E ex PM)/W Mex P(-A) deta[D(l7)] 0(U) F(U) (47) 

where deti[Z?(C/)] denotes the fermion determinant with the zero modes omit- 
ted, and F(U) the remnants of the observable involving fermion fields after 
incorporating the zero modes from the fermion determinant, which in general 
can be expressed in terms of the eigenfunctions of D. It is evident that the rel- 
evant topological sectors which have nonzero contributions to (|47| ) must come 
in ±i/ pairs ( v ^ ). 

Now consider Monte Carlo simulations of lattice QCD with dynamical 
quarks. Suppose the initial gauge configuration is {U} in the topological sec- 
tor of winding number v. A trial configuration {U'} with winding number v' 
is generated. We intend to reject {U'} if u' ^ v. Then, for v' = u, {U'} is 
accepted according to the Metropolis algorithm. That is, we compute the ratio 

P({U}^{U'}) = C ^^eM-MU')+A g (U)) (48) 

and compare P with a random number x drawn from the uniform distribution 
in (0, 1). If x < P, then {U'} is accepted, otherwise {U'} is rejected. Now our 
problem is how to reject {U'} with v' ^ u, before the Metropolis sampling. 

11 



In general, for any two configurations {U} and {U'}, (£48]) can be written 
as 

P({U} - {IT}) = ° W/(I ^ IHH) d^^)] ™P(-MU') + A g (U)) (49) 

where deti denotes the fermion determinant with zero modes omittedf]. It is 
evident that P is zero if \is'\ > \is\, but oo if \is'\ < \is\. So, we can easily 
reject {U'} in these two cases. Next we determine whether is' = v or v' — —v. 
This can be done in the following. In reality, the quarks may not be exactly 
massless. However, if there exists a very light quark ( e.g. u or d quark ) of 
almost zero mass which enters the action in the form mipRipi + m*ipL'4'R ( m* 
is complex conjugate of m ), then the factor Q n n\ v \~\ v \i in Eq. fl4"9|) is replaced 
by 

v ; m v 0{v) + m*- v 6(-v) V ; 

where 6(x) is the step function. In order to reject the configurations {U'} with 
winding number v' ^ ±z/, we first set m* — m ~ 0, and compute P according 
to its definition (f4"8|). Then we reject {U'} if P < Pi(m) or P > p 2 (m), where 
Pi(m) and P2{i r n) are easily determined since P becomes either zero or infinity 
for v' 7^ ±u, in the limit m — > 0. Now we have {U'} with v' = ±v. Next we 
restore m* ^ m and compute P again using fl4"5|). If P is real, then v' = is, 
according to (|50D , so {C/'} is further sampled by Metropolis algorithm. On the 
other hand, if P is complex, then is' = —is, but the configuration {U'} can 
be used for the sector of winding number —is rather than being thrown away. 
As mentiond above, in general, the relevant topological sectors for observables 
in QCD must come in ±is pairs. So, our strategy is to perform Monte Carlo 
simulations for both ±is sectors simultaneously such that optimal efficiency 
can be attained. After enough statistics have been accumulated for ±is sectors, 
then we can move on to another pair of relevant topological sectors by starting 
with a prescribed gauge configuration. Finally the expectation value (O), 
(f47|) , can be obtained by averaging measurements over all relevant topological 
sectors with weight factor exp(iis$). 

On the other hand, if D does not possess exact zero modes, then we cannot 
use the above prescription to decompose the link configuartions into distinct 
topological sectors. Furthermore, if D is topologically improper ( or trivial 
), then its fermion determinant must be very different from its counterpart in 
continuum. Even if we introduce a mass term, the fermion determinant of these 
Dirac operators ( e.g., the standard Wilson-Dirac fermion operator[| ) does not 

6 Here we have assumed that D is normal and thus deti(Z?) is positive definite. 

7 We recall that the Wilson-Dirac fermion operator is topologically trivial, which does 
not possess exact zero modes on a finite lattice. 
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tend to the correct value in the chiral limit, thus the effects of dynamical 
fermion cannot be incorporated properly. Therefore the necessary condition 
for D ( or D' — D(l + RD)~ l ) to have the correct fermion determinant is 
that D is topologically proper. Then the locality of D' would in general lead to 
the correct fermion determinant det(-D') on a finite lattice. However, we still 
cannot formulate a criterion ( such as the locality of a topologically proper D' 
) that can gaurantee the correct fermion determinant on a finite lattice. We 
are sure that, if such a criterion exists, it must be more stringent than that 
for the chiral anomaly, since the former involves all eigenvalues of D while the 
latter only the zero modes. 

A viable scheme to reproduce the continuum fermion determinant on the 



lattice is Narayanan and Neuberger's Overlap formalism |?2| which was in- 
spired by Kaplan's Domain- Wall fermion p3| and Frolov-Slavnov's general- 
ized Pauli-Villars regularization [f23j] with an infinite number of Pauli-Villars 
fields. The lattice Dirac operator proposed by Neuberger || fulfils the crucial 
requirements of reproducing correct chiral anomaly and fermion determinant 



on a finite lattice, as reported in ref. [25 



4 Some Examples 

In this section, we illuminate above discussions by computing the axial anomaly 
and the fermion determinant of a topologically proper D which is constructed 
for the purpose of demonstration, and compare with those of the standard 
Wilson-Dirac fermion operator which is topologically trivial. For simplicity, we 
perform these computations on two-dimensional lattices with U(l) background 
gauge field, and compare the numerical results with the exact solution on the 
torus. It is straightforward to extend these numerical computations to four 
dimensional lattices with non-abelian background gauge fields. 

First, we set up the notations for the Wilson-Dirac fermion operator D w 
with Wilson parameter r w > 0, 

D w = m f + - [ 7/ ,(V* + V M ) - r w V* V M ] (51) 

where rrif is the fermion mass, V M and V* are the forward and backward 
difference operators defined as follows. 

Vt4>{x) = U^(x)il)(x + ji) - ip(x) 
V*^(x) = i/>(x) - UUx - jjL)ip{x - /x) 



From Eq. (|T^), the divergence of axial current for the Wilson-Dirac fermion 
operator is 






d^Jfa)) = 2m f tr[^D- 1 {x,x) 

r f trpKV*^)/^ 1 + ^ 1 (V*V,)](x,x)} (52) 
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The last term in (|5^ ) is the axial anomaly which can be written as 

Ms) = 7E( 4 #(^)T 5 ]- trttfax + fi-fU^x)] 

- tr[D-\x + fi, x)^Ul{x)} - Xt[D~ x {x - jl, x)^U^x - /*)] 
-tT[D-\x,x-p)^XJl{x-p)\ }. (53) 

In the classical continuum limit ( a — > ), (|53|) was shown p6|-p9| to agree 
with the Chern-Pontryagin density. On the other hand, the sum of (|5^) over all 
sites on a finite lattice must be zero in the massless limit, as shown in Eq. (J35|), 
since D w does not have exact zero modes. So, in general, the axial anomaly 
fl53"|) of the Wilson-Dirac fermion operator on a finite lattice must be different 
from the Chern-Pontryagin density. Only in the classical continuum limit when 
the lattice spacing is equal to zero and the number of sites is infinite, then the 
exact zero modes could possibly emerge and the correct axial anomaly might 
be recovered. If this does happen, then we would expect that in the massless 
limit, the sum of (^3|) over all sites must undergo a discontinuous transition 
around the point a = 0, for topologically nontrivial gauge background. 

Next, we consider a topologically proper Dirac operator which has exact 
zero modes satisfying the Atiyah-Singer index theorem. Suppose that one has 
constructed a lattice Dirac operator^ as the following 



D = 5 \lDlD w + D W )[ 3\JDlD w - 2D W (54) 



where D w is the Wilson-Dirac operator with a negative mass parameter mo 
and Wilson parameter r w , 

A i , = -m + -MV*+V M )-r w V*V At ] (55) 

For the moment, we fix m = 1 and r w = 1. Since D satisfies the Ginsparg- 



Wilson relation (|T0|) , the anomaly function of D can be easily obtained from 



Eq. (E3). First, we investigate the properties of D [ Eq. (E$) ] in background 



gauge fields with constant field tensor. As we expect, D possesses exact zero 
modes satisfying the Atiyah-Singer index theorem. However, the anomaly 
function is not constant, though the fluctuations are not very large. The 
deviation of the anomaly function of D from the Chern-Pontryagin density is 
defined as 

Sd = t,Z i AMI (56) 



This is a generalized Neuberger-Dirac operator. 
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where the summation runs over all sites, and N s is the total number of sites, 
and A c (x) is the Chern-Pontryagin density which is equal to 

A c {x) = —e^F^ix) (57) 

Z7T 

on a torus. For constant field tensor with topological charge Q, it becomes a 
constant 

A c (x) = ^- (58) 

In Table 1, we list the deviation of the anomaly function of D, 8d, for topolog- 
ical charge Q = 1 to Q = 8 on a 12 x 12 lattice. Although the deviation Sp in 
the third columm seems to be tiny in each case, it should be regarded as large 
when comparing with the deviation 8d' of a local D' [ Eq. ( pD| ) ] listed in the 
fourth coluum. In fact, the fluctuations of the axial anomaly of D can be seen 
clearly by plotting A n (x) versus each site of the lattice, as shown in Fig. 1. The 
horizontal line at height 1/72 denotes the constant Chern-Pontryagin density, 
A c {x) for Q = 1 on a 12 x 12 torus. The triangles denote the axial anomaly 
of D, and its fluctuations are due to the non- locality of D. The average of the 
axial anomaly of D, of course, must agree with the Chern-Pontryagin density, 
since its index is equal to Q. On the other hand, the axial anomaly [ Eq. fl5H|) 
] of the standard Wilson-Dirac fermion operator D w [ Eq. ( |51]) ] with rrif = 
and r w = 1, which is plotted as squares in Fig. 1, is very different from the 
Chern-Pontryagin density at each site. This is due to the fact that D w does 
not have exact zero modes, thus the sum of its axial anomaly over all sites is 
zero on any finite lattice. So, the deviation of the axial anomaly of massless 
Wilson-Dirac operator, Sjj w , is equal to one for any Q, as listed in the last 
coluum of Table 1. 

The non-locality of D can be easily seen by plotting D(x,y) as a function 
of \x — y\. In Fig. 2, one of the Dirac components of \D(x, 0)|, say, \Du(x, 0)| 
is plotted versus the distance |x|, on a 12 x 12 lattice with periodic boundary 
conditions, in a constant background gauge field with topological charge Q — 1. 
It is obvious that D(x,y) is nonlocal. 

Although D is topologically proper, the non-locality of D must yield in- 
correct chiral anomaly and fermion determinant. In Table 2, we compute the 
fermion determinant of D for different topological sectors of constant field 



tensor, and compare them with the exact solution on a torus [BC 



deti[£>(Q)] ea;act = iV 



\ 




(59) 



where the zero modes have been omitted and the normalization constant is 



fixed by iV = J j—^ such that deti[D(l)] exact — 1. It is evident that the 
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Q 


A c (x) = 2Q/(L 1 L 2 ) 


S D 


5 D > 


fc» 


1 


1.389 x 10- 2 


6.232 x 10~ 2 


2.687 x 10~ 4 


1.0 


2 


2.778 x 1CT 2 


2.434 x 1CT 2 


2.413 x 10~ 4 


1.0 


3 


4.167 x 10~ 2 


8.026 x 10" 3 


2.036 x 10~ 4 


1.0 


4 


5.556 x 1CT 2 


3.192 x lCT 3 


1.982 x 10~ 4 


1.0 


5 


6.944 x 10 -2 


1.026 x lO" 3 


1.363 x 10~ 4 


1.0 


6 


8.333 x 10~ 2 


5.149 x 10~ 4 


1.184 x 10~ 4 


1.0 


7 


9.722 x 10~ 2 


1.910 x 10~ 4 


5.873 x 10" 6 


1.0 


8 


1.111 x 10" 1 


1.061 x 10~ 4 


3.446 x 10~ 5 


1.0 



Table 1: The chiral anomaly versus the topological charge Q. The exact 
solution is listed in the second column. The deviations of the chiral anomaly 
for D [ Eq. ([54]) ] and D' [ Eq. ( |60|) ] are listed in the third and fourth coluums 
respectively. The last column lists the deviations for the Wilson-Dirac fermion. 



fermion determinant of D ( in the third coluum ) disagrees with the exact 
solution ( in the second coluum ), while those of a local D' [ Eq. fl5"0| ) ] ( in the 
fourth coluum ) agree fairly well with the exact solutions. The fermion deter- 
minants for the Wilson-Dirac operator are listed in the fifth coluum, where all 
eigenvalues of D w are taken into account. Thus its fermion determinant must 
become smaller for larger Q due to the emergence of small real eigenvalues. If 
we discard these small real eigenvalues of D w , which correspond to the exact 
zero modes of D, then the resulting fermion determinant of D w is listed in the 
last coluum of Table 2, but still in disagreement with the exact solution. In 
general, the fermion determinant ratio of D w between any two gauge config- 
urations ( no matter within the same topological sector, or in two different 
sectors ) disagrees with the exact solution, except for the trivial sector with 

Q = 0. 

As long as we have a topologically proper D, we can use the topologi- 
cally invariant transformation (P to obtain a local D' such that the anomaly 
function and the fermion determinant both agree with the exact solutions for 
smooth gauge configurations. For simplicity, we fix R = rl in ([3]). Note that 
there is a large range of r values to render D' local for constant background 
gauge fields. For instance, we can pick r = 0.5, then ([[]) gives 



D' = 10 ^DlD 



D„ 



ii\/ dId 



D„ 



(60) 



It is easy to check that D' is local by plotting \D'(x, y)\ versus \x — y\, as shown 
in Fig. 3. The chiral anomaly and fermion determinant of D' agree very well 
with the exact solutions, as listed in Table 1 and Table 2. 

Now we can introduce local fluctuations to the field tensor through the 
sinusoidal terms in Eqs. (|7^) and (|73|). For Q = 1, h\ = 0.1, h 2 = 0.2, 
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Q 


det x [D(Q)] exact 


detiP(Q)] 


defci[D'(Q)] 


det[X>«,(Q)] 


det' [D W (Q)} 


1 


1.00000 


1.00000 


1.00000 


1.00000 


1.00000 


2 


4.24264 


5.21339 


4.22458 


0.34182 


3.98005 


3 


13.8564 


23.5281 


13.7284 


0.14620 


11.8619 


4 


38.1838 


98.6168 


37.2195 


0.07092 


28.8108 


5 


92.7342 


394.114 


90.3215 


0.03841 


61.1240 


6 


203.647 


1535.83 


198.435 


0.02185 


112.466 


7 


411.296 


5759.04 


391.581 


0.01337 


190.058 


8 


773.221 


21275.2 


726.985 


0.00854 


293.285 



Table 2: The fermion determinant versus the topological charge Q. The exact 
solution for a torus of size 12 x 12 is listed in the second column, which is 
computed according to Eq. (|59). The normalization constants are fixed by 
the first row in each case. The lattice Dirac operators D, D' and D w are 
defined in Eqs. (|54D, fl60"D and (|5T|) respectively. 



|(0) 



A\ = 0.3, A 2 = 0.4 and n\ = n 2 = 1 on a 12 x 12 lattice, the anomaly 
function of D' is computed and compared with the Chern-Pontryagin density 
on the torus, as shown in Fig. |j. It is evident that A n (x) agress with A c (x) 
very well at each site. 

Next we consider another lattice Dirac operator D" which has the same 
form of D' ( |5^) except tuq = 0.1 in the D w , 



D" = 10 yJDlD w + D W ) 11 J DlD 



D„ 



(61) 
(62) 



where 



D„ 



-0.l + i[ 7|l (V* + V / .) 



v*,v,] 



Again we investigate the properties of D" in a background gauge field of con- 
stant field tensor. Then we find that D" does not have any zero modes for 
Q > 5, in contrast to the behaviors of D and D' having zero modes satisfy- 
ing the Atiyah-Singer index theorem. In terms of the general index theorem, 
Eq. Q2"T|), we list the topological characteristics of D, D' and D" in Table |^. It 
shows clearly that c[D n ] depends on the gauge configuration through D" . Note 
that D is nonlocal for Q — 1 to Q — 8, but the Atiyah-Singer index theorem 
is still satisfied. We recall that the nonlocal D and the local D' are related by 
the topological invariant transformation (fjj which preserves the zero modes as 
well as the index. 

Now we introduce local fluctuations to the field tensor through the sinu- 
soidal term in Eq. ([721). Starting with constant field tensor at Q — 4, we 
gradually increase the amplitude of the sinusoidal term, A\ with n 2 = 1. The 
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Q 


c[D] 


c[D'} 


c[D"} 


1 


1 


1 


1 


2 


1 


1 


1 


3 


1 


1 


1 


4 


1 


1 


1 


5 


1 


1 





6 


1 


1 





7 


1 


1 





8 


1 


1 






Table 3: The topological characteristics of D, D' and D" versus the topological 
charge Q on a 12 x 12 lattice. 



results are shown in Table [|. For A\ < 0.0802, the index of D" is 4, satisfying 
the Atiyah-Singer index theorem. However, around A\ = 0.0803, the index 
of D" becomes 3, thus the index theorem (|27]) is violated for any integer c[D"]. 
Therefore, we deduce that, in general, if c[D] is not a constant for any integer 
topological charge Q, we can generalize c[D] to a rational number, then the 
index theorem (|27D also holds for any rough gauge fields. It is remarkable that 
the index of D" changes so abruptly around A\ = 0.0802 ~ 0.0803, a singular 
behavior characterizes the topological phase transition. It is very unimagin- 
able for us to see odd number of zero modes ( iV_ = 3 ) in a background gauge 
field of even number of topological charges ( Q — 4 ). The underlying non- 
perturbative mechanism is very intriguing, and might have many unexpected 
consequences. 

We note in passing that if we substitute D' [ Eq. ( |60|) ] and D" [ Eq. fl6l"D 
] by the Neuberger-Dirac operator D h = I + D w (Dl ] D w )^ 1 ^ 2 with m = 1 
and m = 0.1 respectively. The essential features of above numerical results 
are almost the same. In particular, the good agreement between the anomaly 
function and the Chern-Pontryagin density at each site, as shown in Fig. 4. 
We also note that the geometrical aspects of chiral anomalies in the overlap 



was investigated by Neuberger [ |31|j . 



5 Summary and Discussions 

In summary, we have shown that the condition for lattice Dirac operator D 
( or D' = D(l + RD)^ 1 ) to reproduce the continuum anomaly function 
at each site of a finite lattice for smooth background gauge fields is that D 
possesses exact zero modes satisfying the Atiyah-Singer index theorem ( i.e., D 
is topologically proper ) . This is the basic requirement for any Dirac operator to 
be used for lattice QCD if one wishes to reproduce correct chiral anomaly and 
fermion determinant ( ratio ) on a finite lattice. It is evident that the standard 



18 



Q 


4°> 


N_ -N+ 


c[D»] 


4 


0.0100 


4 


1 


4 


0.0200 


4 


1 


4 


0.0400 


4 


1 


4 


0.0600 


4 


1 


4 


0.0800 


4 


1 


4 


0.0801 


4 


1 


4 


0.0802 


4 


1 


4 


0.0803 


3 


3/4 


4 


0.0804 


3 


3/4 


4 


0.0810 


3 


3/4 


4 


0.0820 


3 


3/4 


4 


0.0900 


3 


3/4 



Table 4: The topological characteristics of D" versus the amplitude of the 
local sinusoidal fluctuation on a 12 x 12 lattice, where c[D"] is obtained in 
accordance with the general index theorem 



Wilson-Dirac fermion operator does not meet this basic requirement. If D is 
topologically proper but nonlocal, then one can try to use ([I]) to transform D 
into a local D' = D(l + RD)~ l such that the anomaly function of D' agrees 
with the Chern-Pontryagin density in continuum. Note that there is no need 
for fine-tunings of R. The set of transformations {T(R)} defined in (|XJ) form 
an abelian group with the group parameter space {R}. It is evident that ([!]) 
is one of the simplest topologically invariant transformations which preserve 
the zero modes and the index of D. 

The chiral anomaly on the lattice can be understood from two seemingly dif- 
ferent viewpoints. In general, if we have lattice Dirac operator D which breaks 
the usual chiral symmetry according to Eq. (|TTD, then the chiral anomaly can 
be obtained through the chiral symmetry breaking term [ Eq. (^ ] , while the 
fermion measure is invariant under the usual chiral transformation, in contrast 
to the continuum theory in which the non-invariance of the fermion measure 
under the chiral transformation leads to the emergence of chiral anomaly after 
regularization JJIJ. Another viewpoint to understand the chiral anomaly on a 
finite lattice is only for those D having chirally symmetric limit D c under the 
transformation (|l|). Then D satisfies the Ginsparg- Wilson relation ([H]) which 
is the lattice chiral symmtry, thus the lattice action is invariant under the 
lattice chiral transformation. However, the fermionic measure is not invariant 
under the lattice chiral transformation, and thus leads to the chiral anomaly 
[F21J, reminiscent of the situation in continuum. Therefore, for Ginsparg- Wilson 
Dirac operator, one must obtain the same anomaly function (^) no matter 
which viewpoint one prefers. However, Eqs. Q25| ) and ( |5T] ) are the general 
formulas of the anomaly function for any lattice Dirac operator. 
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There is a vital difference between the continuum Dirac operator and the 
lattice ones. In continuum, the Dirac operator in a smooth non-trivial gauge 
background must have exact zero modes satisfying the Atiyah-Singer index 
theorem. On the other hand, lattice Dirac operator does not necessarily possess 
exact zero modes in a smooth non-trivial gauge background, not to mention 
satisfying the Atiyah-Singer index theorem. This naturally leads to the concept 
of topological characteristics, c[D], which is an intrinsic attribute of D, in 
general, is not entirely due to the non-locality of D and/or the occurrence of 
species doubling. Even if one starts with a lattice Dirac operator D which is 
local and free of species doubling in a constant background gauge field, and the 
zero modes of D also satisfies the Atiyah-Singer index theorem. However, as 
we increase the topological charge or the local fluctuations of the background 
gauge field, D may undergo a topological phase transition and its zero modes 
may violate the Atiyah-Singer index theorem. In fact D has become nonlocal 
long before it reaches the topological phase transition point. So, the non- 
locality of D cannot be the crucial factor for violations of the Atiyah-Singer 
index theorem, though it explains the breakdown of the anomaly function fl29|) 
unambiguously. If one tries to explain the violations of Atiyah-Singer index 
theorem in terms of species doubling, then one has great difficulties to explain 
the occurrence of odd number of zero modes in a gauge background of even 
number of topological charges, as shown in Table f|. Therefore we conclude that 
the topological characteristics of lattice Dirac operator D is a basic attribute 
of D, which in general is a rational number appearing in the general index 
theorem (|27| ) for integer topological charge Q. 

Given a lattice Dirac operator D, we cannot gaurantee that for any gauge 
configuration, there exists an operator R such that the topologically invariant 
transformation (H) gives a local and smooth D' . The notion of locality is 
usually defined as follows : \D(x,y)\ ~ exp(— M\x — y\/a) with M = 0(1), or 
\D(x, y)\ = for \x— y\ > d with d = 0(a). Although \D(x, y)\ can be fitted by 
exp(— M\x — y\/a) with M = 0(1), it may turn out that \D(x, y)\ is still highly 
fluctuating and/or anisotropic. So, it is necessary to introduce the notion of 
"smoothness" which measures the fluctuations of \D(x,y)\ with respect to 
exp(— M\x — y\/a). To determine the criterion for the gauge configuration 
such that D' is local and smooth is beyond the scope of the present paper. 
It is apparent that such criterion should take into account of the fluctuations 
between the plaquettes as well as the deviations of each plaquette from the 
identity. Otherwise the criterion is only useful in the weak field limit. For 
example, the Neuberger-Dirac operator Dh with itiq = r w = 1 and D' [ Eq. 
fl60"|) ] are both local and smooth for gauge configurations with constant field 
tensors up to very large Q values ( see Table 3 ), though each plaquette is 
very different from the identity. If one takes into account of the degree of 
freedom provided by the transformation ([!]), then the criterion for a smooth 
gauge configuration could be less stringent than without it. 
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A scheme for dynamical fermion simulation of lattice QCD is outlined in 
Section 3. Since the direct computation of fermion determinant is prohibitively 
expensive for lattice QCD, a practical implementation of the scheme must be 
worked out before we can proceed to any realistic calculations. Nevertheless, 
the prescription of decomposing link configurations into distinct topological 
sectors and performing Monte Carlo simulations for each relevant topological 
sector individually may open new avenues to tackle the longstanding problems 
in QCD. 



In this appendix, we give an example of prescribed background gauge field 
which has been used in ref. [[25|, [12| for two dimensional and four dimensional 
lattices respectively. On a 4-dimensional torus ( x^ G [0, L J, u — 1, • ■ ■ , 4 



the simplest nontrivial gauge fields can be represented as 



-7*1, 



Ai(ar) = 
A 2 {x) ~- 
A 3 (x) -- 
A 4 (x) -- 


= t a 
= t a 
= t a 
= t a 


"2tt/h 27rq lX2 | 4(0) „./ 
L\ L\Ij2 * 

2nh2 1 l (0) "inf 27rni o;V 


'2im 2 \" 
, L 2 X V\ 

'2im A \" 


~2nh 3 27rq 2 x A ^ 4(0)^/ 
L 3 L 3 L A * 
" 27r ^4 (o) /27rn 3 \ 


, u x v\ 


[ U 4 V L 3 Xa )_ 





(63) 
(64) 
(65) 
(66) 



where t a is any one of the generators of the gauge group with the normalization 
tr(t a t b ) = n 5 ab , and q% and q 2 are integers. The global part is characterized 
by the topological charge 



Q = ^TTT / <^ x e M^<x tr (F^Fxff) = n q x q 2 

6Z7T J 



(67) 



which is an integer. The harmonic parts are parameterized by four real con- 
stants hi, h 2 , h 3 and h A . The local parts are chosen to be sinusoidal fluctuations 



(0) *(0) 4(0) 



|(0) 



A 3 V) and A±' , and frequencies ^f^ 2 , 



2-wn-i 



2-p* and 



with amplitudes A\ , A 2 . . . : . ,,.....,.,, ...... .,,.,,, s L 

^j^- where n\, n 2 , n 3 and n 4 are integers. The discontinuity of A\(x) ( A 3 (x) 
) at x 2 = L 2 ( X4 = L 4 ) due to the global part only amounts to a gauge 
transformation. The field tensors Fi 2 and F34 are continuous on the torus, 
while other F's are zero. 

To transcribe the continuum gauge fields to the lattice, we take the lattice 

iV> 



sites at x M = 0, a, ..., (iV M — l)a, where a is the lattice spacing and L h 
is the lattice size. Then the link variables are 



U\{x) = exp [L4i(x)a] 



(68) 
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U 2 (x) 
U 3 (x) 
UAx) 



iA 2 {x)a + Wxa,(j>r 3 -i)o — = t 



cxp 

exp [L4 3 (o;)a 

exp 



Lx 



2irq 2 x 
iA 4 {x)a + i<3 X4 ,(jv 4 -i)o — j t 



(69) 
(70) 
(71) 



The last term in the exponent of U 2 (x) ( U±(x) ) is included to ensure that 
the field tensor F\ 2 ( F 34 ) which is defined by the ordered product of link 
variables around a plaquette is continuous on the torus. For U(l) gauge field 
on a two-dimensional lattice, the link variabes are reduced to Eqs. (p5|) and 
191) with 



A x {x) 

A 2 (x) 



27r/ii 2tcQx 2 (o) . {2irn 2 

+ A{ sm [ — - — x 2 



2irh 2 



L\L 2 

( o) . /2vmi 
+ A\ sm — - — X\ 
V Lx 



V L 5 



where the topological charge is 



(72) 
(73) 



Q = — I d z x F 12 



(74) 
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Figure 1: The anomaly function on a 12 x 12 lattice. The background gauge 
field has constant field tensor with topological charge Q — 1. The horizontal 
line at height 1/72 is the Chern-Pontryagin density on the torus. The axial 
anomaly of D [ Eq. (]53j) ] is denoted by triangles which fluctuate around the 
exact solution, while that of the Wilson-Dirac operator [ Eq. (|5l| ) with m/ = 
and r w — 1 ] by squares which are far away from the exact solution. The axial 
anomaly of D' [ Eq. fl60|) ] coincides with the exact solution ( the horizontal 
line at height 1/72 ). 
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Figure 2: One of the Dirac components of D(x, 0), \Du(x, 0)|, is plotted as 
a function of |x|. The lattice is 12 x 12 with periodic boundary conditions. 
The background gauge field has constant field tensor with topological charge 
Q — 1. All data points with the same distance \x\ from the origin have been 
averaged. The solid line is an exponential fit to the data points. 
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Figure 3: \D n (x, 0)| is plotted as a function of \x\. The lattice is 12 x 12 
with periodic boundary conditions. The constant background gauge field has 
topological charge Q — 1. All data points at the same distance \x\ from 
the origin have been averaged. The solid line is an exponential fit to the 
data points. The same decay constant also fits very well for all other Dirac 
components of D(x, y) and for any reference point y. 
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Figure 4: The anomaly function of D' [ Eq. ( |60|) ] on a 12 x 12 lattice versus the 
Chern-Pontryagin density A c {x) [ Eq. (|57|) ] on the torus. The parameters of 
the background gauge field are Q = 1, hi — 0.1, h 2 = 0.2, A\ = 0.3, A" 2 = 0.4 
and 7ii = n2 = l- The Chern-Pontryagin density is denoted by circles which 
are connected by line segments for the visual purpose. The anomaly function 
of D' is denoted by triangles. 
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